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Abstract 

Imaging air Cherenkov telescopes (lACTs) detect the Cherenkov hght from exten- 
sive air showers (EAS) initiated by very high energy (VHE) 7-rays impinging on 
the Earth's atmosphere. Due to the overwhelming background from hadron induced 
EAS, the discrimination of the rare 7-like events is vital. The influence of the geo- 
magnetic field (GF) on the development of EAS can further complicate the imaging 
air Cherenkov technique. The amount and the angular distribution of Cherenkov 
light from EAS can be obtained by means of Monte Carlo (MC) simulations. Here 
we present the results from dedicated MC studies of GF effects on images from 7-ray 
initiated EAS for the MAGIC telescope site, where the GF strength is ~ 40 /iT. The 
results from the MC studies suggest that GF effects degrade not only measurements 
of very low energy 7-rays below ~ 100 GeV but also those at TeV-energies. 
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1 1 Introduction 



2 Imaging air Cherenkov telescopes (lACTs) aim at the detection of Cherenkov 

3 light from extensive air showers (EAS) initiated by very high energy (VHE) 7- 

4 rays impinging on the Earth's atmosphere. lACTs make use of the differences 

5 between the angular distributions of Cherenkov light from 7-ray and hadron 

6 induced EAS to efficiently discriminate the hadronic background. Due to the 

7 overwhelming abundance of hadrons in the cosmic rays (mostly protons), the 

8 discrimination of the rare 7-ray events is rather difficult, in particular for 

9 energies below 100 GeV. Only a small fraction of the recorded data are due 

10 to 7-ray initiated EAS. The influence of the geomagnetic field (GF) on the 

11 development of EAS can further complicate the suppression of the hadronic 

12 background and therefore reduces the sensitivity of an instrument. In addi- 

13 tion, GF effects can systematically affect the energy resolution of an lACT 

14 (USE]. 

15 The influence of the GF on EAS was already qualitatively discussed in 1953 0]. 

16 It was pointed out that the east-west separation of electrons and positrons in 

17 EAS due to the GF can be non-negligible compared to the displacement from 
IB multiple Coulomb scattering. Furthermore, the Lorentz force systematically 

19 deflects the particles into opposite directions whereas the displacement due to 

20 multiple Coulomb scattering is random. It was also argued that GF effects are 

21 relatively less important for hadron induced EAS than for 7-ray induced EAS. 

22 The scattering angles occurring in nuclear interactions of hadronic EAS give 

23 rise to a lateral displacement of the shower particles much larger than that 

24 due to the influence of the GF. The comparatively large transverse momenta 

25 of the secondary particles result in a large angular spread of the directions of 

26 the electromagnetic sub-cascades generated by pion decay. 

27 The influence of the GF on the average lateral spread of atmospheric Cherenkov 

28 radiation was studied by means of computer simulations already more than 

29 30 years ago [5j . GF effects on real 7-ray and proton initiated EAS were later 

30 on studied using a non-imaging Cherenkov telescope [6]. It was reported that 

31 the influence of the GF on proton initiated EAS results in a significant reduc- 

32 tion of the count rate. It was shown elsewhere [7J that lACT measurements 

33 of TeV 7-rays from the Crab nebula were not significantly affected when the 

34 component of the GF normal to the shower axis, i.e. transversal component of 

35 the GF, was below 35 /iT. The average shape and reconstructed intensity of 

36 Cherenkov images from hadrons was found to be independent of the transver- 

37 sal component of the GF. However, it was pointed out that the instrument 

38 was not sensitive enough to study GF effects. More recent measurements of 

39 7-ray showers carried out with a transversal component of the GF strength of 

40 \B±\ > 40 /iT revealed GF effects in observational data compatible to those 

41 predicted by MC simulations, both for 7-ray and hadron showers p[9] . The 

42 authors suggest that for EAS developing under unfavourable orientation with 

43 respect to the direction of the GF the corresponding Cherenkov light images 
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44 in the camera of an lACT will be rotated. As the information on the ori- 

45 entation of shower images provides the most powerful discrimination between 

46 7-ray shower images from a point-like source and any unwanted isotropic back- 

47 ground (mainly due to hadrons) this results in a degradation of the sensitivity 

48 of an lACT. However, it was demonstrated that a correction for GF effects in 

49 7-ray initiated Cherenkov images is possibly resulting in an increased detec- 

50 tion significance and better sensitivity of the lACT [T OfTT] . The correction for 

51 GF effects required simulated 7-ray showers. 

52 lACTs currently in operation offer improved imaging capabilities, i.e. better 

53 optical point spread function (PSF), pixel resolution and timing capabilities 

54 of the electronics, and are therefore more sensitive to GF effects than previous 

55 instruments. 

56 In this paper we present results from dedicated MC studies of GF effects on the 

57 imaging technique. The studies were carried out for the Major Atmospheric 

58 Gamma-ray Imaging Cherenkov (MAGIC) telescope [T^lIT^ . which is located 

59 on the Canary Island of La Palma at the Roque de los Muchachos Observatory 

60 at 2200m ahitude (28.45° N,17.54° W). 



61 2 The MAGIC Telescope 



62 The 17 m diameter MAGIC telescope is currently the largest single dish lACT 

63 in operation. The imaging camera in the focal plane of the tessellated parabolic 

64 reflector consists of 577 photomultiplier tubes (PMTs), all of which are ar- 

65 ranged in a hexagonal configuration. The inner part of the camera is equipped 

66 with 397 PMTs of a diameter of 0.1° whereas the outer part is equipped with 

67 larger PMTs of a diameter of 0.2°. The reflector has a focal length of 17m. 

68 The field of view of the camera is 3.5° and the angular resolution for 7-rays 

69 is about 0.1°, depending on the energy. The telescope is in continuous oper- 

70 ation since summer 2004. It allows for a detection of a 7-ray source with an 

71 absolute intensity of ~ 2 % of the Crab nebula and similar energy spectrum 

72 within 50 hours at energies > 200 GeV on a significance level of 5 standard 

73 deviations. 

74 MAGIC is currently being upgraded through the addition of a twin telescope 

75 to achieve an improved sensitivity and a lower energy threshold [H]. Further 

76 technical details and information on the performance of the instrument can 

77 be found elsewhere [T5] . 
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78 3 The Geomagnetic Field at the MAGIC Site 



79 The component of the GF normal to the shower axis is relevant for the east- 

Bo west separation of electrons and positrons during the shower development. 

Bi For this study the telescope optical axis has always been set parallel to the 

B2 direction of the primary 7-ray. 




Fig. 1. (a) The absolute value of the component of the GF normal to the direction 
of the EAS versus azimuth (Az) angle and zenith angle (ZA) for the Roque de los 
Muchachos observatory on La Palma. (b) The definition of the coordinate system 
used throughout this work. 9 denotes the angle between the direction of the EAS and 
the direction of the GF. The Az angle is defined like in the CORSIKA program |16j . 
i.e. it refers to the momentum of the incoming 7-ray and is counted counterclockwise 
from the positive x-axis towards west. The telescope optical axis has always been 
set parallel to the direction of the primary 7-ray. 

83 Figure [1] (a) shows the absolute value of the GF component |-B_l| normal to the 

84 direction of the EAS versus azimuth (Az) angle and zenith angle (ZA) for the 

85 MAGIC telescope site. The value was determined for 10 km a.s.l. according 

86 to the epoch 2005 International Geomagnetic Reference Field (IGRF) model 

87 [IT]. The MAGIC telescope is focused to a distance of 10km a.s.l., which is 

88 the most likely location of the shower maximum for 100 GeV 7-ray induced 

89 EAS at small ZAs. The trajectory of the strongest source of steady VHE 7-ray 

90 emission in the Galaxy, the Crab nebula, is indicated. As the magnetic field 

91 lines at La Palma are tilted by ~ 7° westwards with respect to the meridian 

92 [17] the trajectory is asymmetric with respect to 180° Az angle. 

93 For La Palma, the minimum influence of the GF is expected to occur for EAS 

94 developing in direction of the magnetic north at ZA = (90° — I) ^ 51° and 

95 Az = 0°, where the angle 6 between the shower axis and the GF becomes 

96 smallest (see figure [1] (b)), i.e. for EAS developing along the field lines. / 

97 denotes the angle under which the GF lines dip into the Earth's surface, which 
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98 is ~ 39° for La Palma [T7] . Hence, the maximum influence is expected for EAS 

99 developing perpendicular to the direction of the GF lines, i.e. for ZA ^ 39° 

100 and Az = 180°. 

101 The results from the MC studies on the GF effects presented here are speciflc 

102 to the MAGIC telescope because of the local GF strength and the relevant 

103 telescope parameters like the reflector area, camera pixelisation and the 7-ray 

104 PSF. 



105 4 Monte Carlo Simulations 



106 The production of Monte Carlo (MC) data for MAGIC involves three steps 

107 [18] : 

108 1. The CORSIKA (COsmic Ray Simulations for KAscade) MC program (ver- 

109 sion 6.019) [16] is used to simulate the development of 7-ray and hadron 
no induced extensive air showers (EAS) and the production of Cherenkov light 

111 for a given set of input parameters, like the primary 7-ray or hadron en- 

112 ergy, height above sea level, the magnitude and direction of the GF, and 

113 so on. The GF components are set to the values for the location of the 

114 MAGIC telescope (La Palma, 28.8° N, 17.9° W) which are provided by the 

115 IGRF model [T7j. The x-axis of the Cartesian CORSIKA reference frame is 

116 aligned with the magnetic north pole and the |/-axis points to the west. The 

117 Az angle is counted counterclockwise from the positive x-axis and refers to 

118 the direction of the primary 7-ray (flgured] (b)). Furthermore, the so-called 

119 US standard atmosphere is used as a model for the Earth's atmosphere. 

120 2. The binary output of CORSIKA, containing information on the Cherenkov 

121 photon direction and its position on ground, is processed with a dedicated 

122 Reflector program, which does the ray-tracing of the Cherenkov photons. 

123 To be able to adapt to different conditions without being forced to rerun 

124 CORSIKA, atmospheric absorption and scattering of Cherenkov photons as 

125 well as mirror condition is taken into account at this stage. 

126 3. Finally, the output of the Reflector program is processed by the Camera 

127 program simulating the entire readout chain, i.e. PMT response, trigger 

128 and data acquisition system including electronic noise. Normally, a compact 

129 next-neighbour coincidence is required, i.e. at least four neighbouring pixels 

130 are required to trigger, and, if any of the pixels is taken out of the group, 

131 the remaining pixels are still neighbours. To adapt the MC data to the 

132 optical performance of the telescope, the simulation of the optical point 

133 spread function (PSF) can be tuned at this stage. The calibration and the 

134 image parameter calculation (Hillas analysis [19]) is done using the MAGIC 

135 Analysis and Reconstruction Software (MARS) [20] . 
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136 4-1 Dedicated MC Production 

137 For the present studies only 7-rays were simulated. The MC data were pro- 

138 duced following for most instances the standard MC production of the MAGIC 

139 telescope as described beforehand. All events were simulated as originating 

140 from a point source. By definition, the telescope optical axis is always parallel 

141 to the direction of the primary 7-ray. The impact parameter (IP) is defined as 

142 the distance from the centre of the telescope mirror to the shower axis, which 

143 has the same direction as the primary 7-ray. 



y (West) 




(a) (b) 

Fig. 2. (a) An illustration of the setup used in the MC simulations: the telescope 
positions are indicated as full circles and the FAS core location (the primary 7-ray 
impact point on the ground) as an open circle. The impact distance on ground r is 
defined as the distance between the FAS core location and the telescope position. 
The telescope is situated always somewhere in the circle within which Cherenkov 
photons are kept. The simulations were done for fixed impact distances on ground 
r = 20 m,. . . 180 m and angles (p = 0°, 30°,. . .330°. (b) The relation between the 
impact parameter (IP) and the impact distance r on the ground is illustrated. In 
direction of the telescope's inclination the impact parameter equals rcos(ZA). 

144 In contrast to the production of standard MC data, where the EAS core loca- 

145 tion is randomly placed somewhere in a circle on the plane perpendicular to the 

146 direction of the EAS (to estimate the telescope effective collection area), the 

147 EAS for this study were simulated for fixed core locations and all Cherenkov 

148 photons arriving in a circle of 200 m radius were kept. This procedure reduces 

149 computing time because each CORSIKA event can be used multiple times by 

150 placing the telescope (at the level of the Reflector program) somewhere into 

151 the circle within which Cherenkov photons are kept (figure [2] (a)). Besides, 

152 this approach allows to study the influence of the OF on the shower images 

153 in great detail. For technical reasons the telescope was placed on equidistant 

154 points concentrically to the EAS core location. By the choice of this setup, 

155 the impact parameter varies like r-i/cos^(Az — (/?)(cos^(ZA) — 1) -|- 1 (figure [2] 
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156 (b)), where r is the distance between the EAS core location and the telescope 

157 position. 

158 The energy of the primary 7-ray was for different samples set to 30 GeV, 

159 50 GeV, 70 GeV, 120 GeV, 170 GeV, 300 GeV, 450 GeV or 1 TeV, respectively. 

160 The ZA was varied between 0° and 60° in steps of 20°, and the Az angle 

161 between 0° and 180° in steps of 30°, as the absolute value of the GF com- 

162 ponent normal to the EAS direction is symmetric in the Az angle (figure [1] 

163 (a)). Hence, the maximum value for the angle 9 achieved with the simulated 

164 telescope orientations is 87°. The choice of discrete values for the 7-ray energy 

165 allows to investigate the energy dispersion of the showers due to the GF. The 

166 distance r between the telescope position and the EAS core location (impact 

167 point of the primary 7-ray on the ground) was varied between 20 m and 180 m 

168 in steps of 20 m and the angle ip (figure [2] (a)) between 0° and 330° in steps 

169 of 30°, resulting in 108 configurations. 

170 About 10^ events were simulated for each 7-ray energy, ZA and Az angle. As 

171 a reference, MC data were also produced without GF. To be as realistic as 

172 possible the MC simulations include the effects of photons from the diffuse 

173 night sky background of 1.75 ■ 10^^ phm~^ s~^ sr~^ at the MAGIC site [21] as 

174 well as electronic noise. 



175 5 Image Analysis 

176 The MC-generated 7-ray showers were analysed using the standard MAGIC 

177 software MARS [20]. Before parameterisation of the shower images a tail-cut 

178 image cleaning was applied (figure [3] (a)) [T5]. The image cleaning requires 

179 the signals to be above a certain level. For the MC studies presented here the 

180 minimum required pixel content was 7photoelectrons (phe) for so-called core 

181 pixels and 4 phe for boundary pixels. 

182 Shower images from 7-ray showers processed with the image cleaning are nar- 

183 row and point towards the source position in the field of view. To a first 

184 approximation the shower images are elliptical and can be described by so- 

185 called Hillas parameters [T^j. Detailed reviews on the imaging technique can 

186 be found elsewhere [7112^ . Some of the image parameters used for these studies 

187 are illustrated in figure [3] (b). The parameters WIDTH and LENGTH char- 

188 acterise the lateral and longitudinal spread of the shower images (minor and 

189 major axes of the so-called Hillas ellipse). Both parameters are very important 

190 since they allow a powerful discrimination of 7-ray images against hadron in- 

191 duced images. The parameters DIST and ALPHA are related to the position 

192 and orientation of shower images in the camera. The parameter DIST is di- 

193 rectly related to the impact parameter of the primary 7-ray (figure [2] (b)). The 

194 image parameter ALPHA is commonly used by standalone lACTs to extract 

195 the 7-ray signal. ALPHA denotes the angle between the major axis of the 
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196 shower image and the vector connecting its centre of gravity with the source 

197 position in the camera plane (camera centre). It provides a very powerful dis- 

198 crimination between 7-ray images from a point-like source and any isotropic 

199 background (mainly due to hadrons), i.e. orientation discrimination. The 7- 

200 ray signal from a VHE 7-ray source under study appears as an excess at small 

201 values in the ALPHA parameter distribution. 




(a) (b) 

Fig. 3. (a) A MC-generated 7-ray shower image after application of the image clean- 
ing. Only pixels surviving the image cleaning are used for the following steps of the 
analysis, (b) Definition of the image parameters. The light distribution is described 
by an ellipse whose major and minor axes describe the longitudinal and lateral 
spread of the Cherenkov light distribution. 

202 Another viable image parameter is the so-called SIZE, which corresponds to 

203 the total integrated light of a shower image after treatment with the image 

204 cleaning procedure. It is therefore an estimate for the primary 7-ray energy. 

205 The so-called DISP method [251121] allows to reconstruct the arrival direction 

206 of 7-ray candidates making use of the shape of a shower image. The DISP 

207 parameter is determined according to the formula 



WIDTH 

DISP = ci(SIZE) + C2(SIZE) 



LENGTH 



208 Therein ci and C2 are second-order polynomials optimised on MC simulated 

209 7-ray showers [25] : 

ci = +1.163° + 0.542° (logio(SIZE) - 2) - 0.672° (logio(SIZE) - 2f , 
C2 = -0.265° - 2.905° (logio(SIZE) - 2) + 2.220° (logio(SIZE) - 2f . 



210 In case of a single telescope the DISP method provides two possible solutions 

211 for the source position. To overcome this ambiguity the asymmetry of the 
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212 shower image along the major image axis, which is related to the longitudinal 

213 development of an EAS in the atmosphere, is used to reconstruct the true 

214 source position. However, owing to false head-tail assignment the percentage 

215 of correctly reconstructed events is typically limited to ~ 80 %, depending on 

216 the 7-ray energy [25]. The outcome of the DISP analysis is usually displayed 

217 in terms of a sky map of arrival directions. Because the DISP parameter 

218 depends on both the eccentricity WIDTH/LENGTH and on the orientation 

219 of the shower images the influence of the GF on the shower development is 

220 expected to degrade also the DISP method. 

221 In this paper, we restrict ourselves to selected but representative results from 

222 the MC study. 



223 6 Results and Discussion 

224 6.1 GF Effects on Shape and Orientation of'j-ray Shower Images 

■ > X (North) 




(a) Az = 0°,e = 12°. (b) Az = 180°, 9 = 87°. 



Fig. 4. Shower images in the telescope camera for 450 GeV 7-rays, 100 m impact 
parameter, ZA = 40°, Az = 0° and 180°, respectively. The images indicated by 
solid lines were obtained for enabled GF and the ones indicated by dashed lines 
were obtained for disabled GF. The orientation of the GF component normal to the 
direction of the EAS (telescope pointing direction) is indicated in the lower right 
part of the figures (see text for more details) . 

225 The orientation of the 7-ray images from a source under study is used for the 

226 suppression of the isotropic hadronic background. It is accordingly important 

227 to study how the GF influences the orientation of shower images. Figure H] 

228 shows shower images (Hillas ellipses) in the telescope camera for 450 GeV 7- 

229 rays, 100 m impact parameter, ZA = 40°, and different Az angles between 0° 
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and 180°. The ellipses drawn with solid lines were obtained for enabled GF 
and the ones drawn with dashed lines for disabled GF in the MC simulation. 
For each angle <y9 = 0°, 30°,. . . 330° the size, the position and the orientation of 
the ellipse in the camera was determined by taking the mean values from the 
corresponding Hillas parameter distribution. The ellipses are superimposed 
on the projection of the CORSIKA coordinate system, whose x-axis is aligned 
with the magnetic north. The direction of the GF component normal to the 
direction of the EAS is indicated in the lower right part of the figures. 
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Fig. 5. The rotation angle of 7-ray images versus energy for 40 m and 120 m im- 
pact parameter, ZA = 0°, 40°, Az = 0° and 180°. Full triangle down data points 
correspond to the case where the connecting line between the EAS core location on 
ground and the telescope position is parallel to the north-south direction (x-axis), 
while for full square data points the telescope is situated on the y-axis. Full triangle 
up and full circle data points correspond to intermediate telescope positions. 
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Fig. 6. The rotation angle of 7-ray images versus impact parameter for ZA = 20° 
and 40°, Az = 0° and 180°, 120 GeV, 450 GeV and 1 TeV energy (see text for more 
details) . 
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238 As can be seen the average orientation is preserved only for shower images 

239 oriented either parallel or normal to the projected direction of the GF in the 

240 camera, which is by definition the GF component normal to the direction of 

241 the EAS (telescope pointing direction). Shower images situated at intermedi- 

242 ate angles are systematically rotated away from the projected direction of the 

243 GF. The sideways spread of the images result in a systematic rotation away 

244 from the camera centre (source position). This was also reported in [TO]. The 

245 magnitude of the rotation depends not only on the angle 6 between the axis 

246 of the EAS and the direction of the GF but mainly on the core position of the 

247 EAS with respect to the telescope, which is given by the angle (p. 

248 By comparing the orientations for images generated with disabled GF to the 

249 ones for enabled GF it is possible to determine the rotation angle. Figure [5] 

250 shows the rotation angle of 7-ray images versus energy for an impact param- 

251 eter of 40 m and 120 m. The figures show that the rotation angle depends on 

252 the impact parameter and on the 7-ray energy. However, the average rotation 
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253 angle for images oriented either parallel or normal to the direction of the GF in 

254 the camera is zero (full square and full triangle down data points). Images ori- 

255 ented at intermediate angles are systematically rotated. For ZA = 0°, Az = 0° 

256 and 40m impact parameter (figured (a)) and for ZA = 40°, Az = 180° (figure 

257 [5] (c) and (d)), which is the most unfavourable telescope pointing direction, 

258 the rotation angle is maximal for 7-ray energies around 450 GeV. 

259 Figure [6] illustrates the dependency of the rotation angle on the impact pa- 

260 rameter. Small impact parameters correspond to images with low eccentricity 

261 WIDTH/LENGTH, which may be rotated through a large angle. Figure [H] (b) 

262 also shows that the direction of the rotation depends on the impact parame- 

263 ter. The shower images are not always rotated away from the direction of the 

264 GF in the camera but can even be rotated towards it. Consequently, the cor- 

265 rection of observational data for GF effects by de-rotating the shower images 

266 must take into account the energy dependence of the rotation angle and its 

267 dependence on the impact parameter (DIST). For the most unfavourable tele- 

268 scope pointing direction and for small impact parameters even 1 TeV shower 

269 images are rotated through a large angle (figure [6] (f)). The influence of the 

270 GF on the shape of the 7-ray shower images also depends on the primary 

271 7-ray energy, the impact parameter and the orientation of the EAS relative to 

272 the direction of the GF. Figure [7] shows the average WIDTH and LENGTH 

273 of 7-ray images versus energy for 120 m impact parameter, ZA = 40°, and 0° 

274 as well as Az = 180°. For the most unfavourable telescope pointing direction, 

275 i.e. ZA = 40° and Az = 180° {9 = 87°) significant GF effects on the image pa- 

276 rameter WIDTH occur for 7-ray energies above ~ 100 GeV (figure [7] (b)). The 

277 influence of the GF on the image parameter WIDTH is considerably larger 

278 than on LENGTH. Due to the influence of the GF the average WIDTH is 

279 increased for images where the connecting line between the EAS core location 

280 on ground and the telescope position is parallel to the magnetic north-south 

281 direction (telescope situated on the x-axis, full and open circle data points). 

282 Images aligned with the direction of the GF in the camera are horizontally 

283 stretched compared to the situation of disabled GF in the MC, whereas im- 

284 ages oriented normal to the direction of the GF have a smaller WIDTH and 

285 are thus elongated due to the influence of the GF (telescope situated on the 

286 y-axis, full and open square data points). 

287 6.2 GF Effects on the Image Parameter ALPHA 

288 In the preceding section it was shown that the GF can strongly alter the 

289 average shape and orientation of 7-ray shower images in the camera. Even 

290 though the orientation of shower images and the image parameter ALPHA 

291 are correlated it is important to investigate the influence of the GF on the 

292 image parameter providing the most powerful discrimination between 7-rays 

293 from a point-like source and unwanted isotropic background (mainly hadrons) . 
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(a) IP = 40 m. (b) IP = 120 m. 

Fig. 8. Normalised distributions of the image parameter ALPHA for 50 GeV 7-rays, 
ZA = 0°, Az = 0° {e = 52°), 40m (a) and 120m impact parameter (b). Different 
configurations are compared: the distributions indicated by sohd hnes correspond 
to if = 90° and 270°, the dash-dotted distributions to tp = 0° and 180°, and the 
distributions indicated by dashed lines to intermediate telescope positions 99 = 30°, 
60°, 120°, 150°, 210°, 240°, 300° and 330°. The corresponding distribution obtained 
without GF in the MC simulation is also plotted (dotted line). The percentage of 
events with |ALPHA| < 9° is given in the legend. 



294 Figures [H]-[in] show the ALPHA distributions (normalised to the number of 

295 entries) for 7-ray energies of 50 GeV, 450 GeV and 1 TeV. Showers recorded at 

296 40 m and 120 m impact parameter were considered together with the different 

297 possible configurations: the distributions indicated by dash-dotted lines cor- 

298 respond to = 0° and 180° (as defined in figure [2] (a)), where the connecting 

299 line between the shower axis and the telescope optical axis is parallel to the 

300 north-south direction. This configuration corresponds to shower images which 

301 are oriented parallel to the direction of the GF in the camera. The distri- 

302 butions indicated by solid lines were obtained for ip = 90° and 270°, where 

303 the connecting line between the shower axis and the telescope optical axis 

304 is parallel to the east-west direction. In this case the shower images are on 

305 average not rotated and oriented normal to the direction of the GF in the 

306 camera. The ALPHA distributions indicated by dashed lines belong to inter- 

307 mediate telescope positions (p = 30°, 60°, 120°, 150°, 210°, 240°, 300° and 

308 330°. The corresponding ALPHA distributions obtained without GF in the 

309 MC simulation are also plotted (dotted lines). The percentage of events with 

310 I ALPHA I < 9° is given in the legend. The cut is indicated by the vertical 

311 dotted line. 

312 It can be seen that for configurations where the connecting line between the 

313 telescope and the shower axis is parallel to the north-south direction (parallel 

314 to the direction of the GF) the corresponding ALPHA distributions (red his- 

315 tograms) can be significantly broadened although the corresponding shower 
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316 images are on average not rotated. However, the ALPHA distribution for the 

317 opposite configuration (connecting fine between the telescope and the shower 

318 axis parallel to the east-west direction) are stronger peaked at low values due 

319 to the infiuence of the GF (compare figure [HI (a) and (b)). The remaining con- 

320 figurations always lead to broadened ALPHA distributions due to the rotation 

321 of the shower images with a preferential direction. 

322 In conclusion it can be stated that the infiuence of the GF can significantly 

323 degrade the orientation discrimination of shower images. It is evident that for 

324 some configurations discussed above the 7-ray signal cannot be recovered by 

325 de-rotating the shower images. 

326 Figure [ID shows the ALPHA distributions for 450 GeV 7-rays, ZA = 40°, 

327 Az = 180° and 40 m as well as 120 m impact parameter. The ALPHA dis- 

328 tributions obtained without GF (dotted line) are shown together with the 

329 distributions obtained for enabled GF (dashed line) and the ones obtained 

330 after de-rotation of the shower images (solid line). It is possible to correct for 

331 GF effects by de-rotating the shower images. However, if all telescope posi- 

332 tions are taken into account (</? = 0° . . . 330°) the improvement in terms of the 

333 percentage of events with |ALPHA| < 9° is less than 8% (figure [TT] (a) and 

334 (b)). Ignoring the most unfavourable configurations with respect to the infiu- 

335 ence of the GF {ip 7^ 0° and 180°) results in ALPHA distributions which are 

336 stronger peaked at small values (figure [TT] (c) and (d)). As expected, for the 

337 most unfavourable configurations {ip = 0° and 180°) the 7-ray signal cannot 

338 be recovered by de-rotation. The corresponding images are not rotated but 

339 their angular distribution is broadened (figure [TT] (e) and (f)). 

340 Given that the energy, the rotation angle and the impact parameter is well 

341 known in MC, the amount of recovered real 7-ray showers from observational 

342 data by de-rotation is expected to be lower. Both the energy and the impact 

343 parameter have to be estimated and are thus known less precisely. To be effi- 

344 cient, the de-rotation of the shower images requires a precise knowledge of the 

345 impact parameter and the information on the energy of the 7-ray candidates 

346 from observational data. Moreover, we focused on intermediate 7-ray energies 

347 where the rotation angle is large (figure [5] (d)) and the spread of the AL- 

348 PHA distribution is rather low. At lower energies than those considered here 

349 the recovery of the 7-ray signal is even less efficient ^J. This is also the case 

350 for 1 TeV 7-rays, where the ALPHA distribution is stronger peaked at small 

351 values and the rotation angle is smaller (figure . Because of the relatively 

352 poor knowledge of the impact parameter in case of real shower images the 

353 improvement in sensitivity by de-rotation of the shower images is expected to 

354 be below 10 %. 
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As figure El but for 450 GeV 7-rays, ZA = 40°, Az = 0° and 180°. 



16 



ISUO, 9 = 90°, 270°, 68.6% 
l§UO, <p = 0M80°, 66.1% 
IbUo, tp?^0°, 90°, 180°. 270°, 
IbUo, allip, 67.3% 



% 10 20 30 40 50 60 70 80 90 

|ALPHA| [°] 



« 0.4 
£0.35 
0.3 
0.25 
0.2 
0.15 
0.1 
0.05 



|B|*0, (p = 90",270°. 83.8% 
l§l*0, (p = 0", 180°, 35.0% 
IbI*0, ((>#0",90", 180", 270", 36.2% 
IbI = 0, all(p, 67.3% 




10 20 30 40 50 60 70 80 90 
|ALPHA| [°] 



Az = 0°, IP 40m, = 12°. 



(b) Az = 180°, IP « 40m, e = 87°. 



IbUo, ip = 90°, 270°, 100.0% 
I SI 7^0, ip = 0°, 180°, 100.0% 
lSl/0, <p7^0°, 90°, 180°, 270°, 100.0% 
ISI = 0, allcp, 100.0% 



% 10 20 30 40 50 60 70 80 90 

|ALPHA| n 



« 0.9 

.0) 

I 0-8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



l§l*0, (p. 90", 270°. 97.1% 
l§l*0, (p = 0°, 180°, 89.1% 
IbI*0, (p*0°,90°, 180°, 270", 90.6%, 
IbI = 0, all*, 100.0% 



t4-|.r:..i.J:r'r 1 



% 10 20 30 40 50 60 70 80 90 

|ALPHA| [°] 



(c) Az = 0°, IP ^ 120m, e = 12°. 



(d) Az = 180°, IP ^ 120 m, 9 = 87°. 



Fig. 10. As figure El but for 1 TeV 7-rays, ZA = 40°, Az = 0° and 180°. 
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Fig. 11. Normalised distributions of the image parameter ALPHA for 450 GeV 
7-rays, ZA = 40° and Az = 180°, 40m ((a), (c) and (e)) and 120 m impact param- 
eter ((b), (d) and (f)). The distributions indicated by dotted hnes were obtained 
without GF and the ones indicated by dashed hnes for enabled GF in the MC 
simulation. For the distributions indicated by solid lines the shower images were 
de-rotated (see text for more details). 
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355 6.3 GF Effects on the DISP -reconstructed 'y-ray Arrival Direction 



356 A DISP analysis of the MC 7-ray shower images was performed to study the 

357 GF effects on the reconstructed arrival directions. 




(a) Az = 0°,e = 12°. (b) Az = 180°, 9 = 87°. 

Fig. 12. DISP-reconstructed arrival directions for 450 GeV 7-rays, impact parame- 
ters between 60 m and 140 m, ZA = 40°, Az = 0° and 180°, respectively (see text 
for more details). 

358 For the optimisation of the second-order polynomials in equation ([T]), dedi- 

359 cated MC 7-ray samples with continuous impact parameter distribution be- 

360 tween m and 500 m were produced. The MC samples were produced for the 

361 same 7-ray energies, ZAs and image cleaning levels as the MC data used for 

362 the preceding studies. The EAS core location was randomly placed somewhere 

363 in a circle on the plane perpendicular to the direction of the EAS. Also, the 

364 MC samples were produced without GF, thus only for Az = 0°. In this way the 

365 results from the DISP method obtained for different telescope pointing direc- 

366 tions are comparable since the DISP polynomials themselves are not subject 

367 to GF effects. 

368 Figure [12] shows the DISP-reconstructed arrival directions for 450 GeV 7-rays, 

369 impact parameters between 60 m and 140 m, ZA = 40°, Az = 0° and 180°. 

370 The projected direction of the GF is indicated in the lower right part of the 

371 figures, and the ellipticity of the distributions of DISP-reconstructed arrival 

372 directions is shown in the lower left part of the figures. The semi-minor and 

373 the semi-major axis of the ellipse correspond to the sigma of a Gaussian fit to 

374 the distributions using bands of Ax,y = ±0.035° parallel and perpendicular 

375 to the projected direction of the GF. The size of the bands was arbitrarily 

376 chosen, but it is not critical for the result. Within errors, the relative differ- 

377 ence of the results obtained from the two orientations is independent of the 

378 size of the bands. The result from the Gaussian fit is also shown in the legend. 

379 As can be seen from the figures the distributions appear to be significantly 
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3B0 elongated perpendicular to the projected direction of the GF, while the peak 
3B1 of the DISP distribution is always centred at the nominal source position 
3B2 (camera centre). The extent of the elongation depends on the angle 9 between 
383 the shower axis and the direction of the GF. The GF effects on the DISP 
3B4 method thus result in a degradation of the sky maps in a way that a point-like 
3B5 7-ray source appears to be extended, i.e. the 7-ray PSF is degraded. The star- 
see shaped appearance of the DISP distributions arises from events with wrong 
3B7 head-tail assignment. The false head-tail assignment cannot be attributed to 
3B8 GF effects since it occurs also for the favourable telescope pointing direction 
3B9 (figure [12] (a)). This illustrates the basic limitation of a single telescope to 

390 properly reconstruct the true source position. In case of a single telescope the 

391 DISP method has to rely on the shower asymmetry along the major axis of 

392 the shower image (see section [5]). 
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Fig. 13. Lateral and longitudinal spread of the DISP distribution versus 7-ray energy 
for impact parameters between 60 m and 140 m, ZA = 40°, Az = 0° and 180°, 
respectively (see text for more details). 



393 Figure [13] shows the lateral and longitudinal spread versus 7-ray energy for 

394 different orientations of the telescope. The spread is defined as the sigma 

395 of a Gaussian fit to the DISP distribution using bands of Ax,y = ±0.035° 

396 parallel and normal to the projected direction of the GF in the camera. Impact 

397 parameters between 60 m and 140 m were considered. The figures clearly show 

398 that, compared to the case of disabled GF, the spread of the DISP distribution 

399 increases significantly for an unfavourable telescope orientation (large angle 

400 9). The maximum spread occurs always perpendicular to the direction of the 

401 GF in the camera. Therefore, depending on the orientation of an EAS with 

402 respect to the telescope and the impact parameter, the DISP-reconstructed 

403 incoming direction of the corresponding primary 7-ray has a large uncertainty. 



20 



404 6-4 GF Effects on the Energy Reconstruction of -ray Images 



405 It was previously discussed that the influence of the GF on the shower de- 

406 velopment affects also the energy reconstruction and the trigger efficiency for 

407 7-rays. 
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(c) ZA = 40°. (d) ZA = 60°. 

Fig. 14. The mean of the image parameter SIZE and the 7 efficiency versus Az angle 
for 450 GeV 7-rays and ZA between 0° and 60°. The telescope is always situated at 
angles (p = Az -|- 90° and 120 m impact parameter. 

408 The east-west separation of electrons and positrons in EAS due to the GF 

409 modifies the Cherenkov distribution on the ground such that the reconstructed 

410 Cherenkov light, i.e. the integrated light content of the shower images is re- 

411 duced for unfavourable shower orientations with respect to the direction of 

412 the GF. 

413 Figure [14] shows the average reconstructed light content of shower images to- 

414 gether with the 7 efficiency versus Az angle for 450 GeV energy 7-rays and 
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415 ZAs between 0° and 180°. To compare equivalent configurations (see figure [2] 

416 (a)) the telescope is always situated at angles ip = kz + 90° and the impact 

417 parameter was set to 120 m. The angle 9 between the direction of the EAS and 

418 the GF is given on top of the abscissa. The 7 efficiency is defined as the ratio 

419 of the number of 7-ray showers surviving the trigger and the image cleaning 

420 to the number of generated 7-rays. From figure [H] (d) it can be seen that the 

421 7 efficiency varies by up to 50 %. For showers close to the trigger threshold 

422 the Cherenkov light distribution on the ground can be thinned out such that 

423 most of the events do not survive the trigger level. This occurs only for show- 

424 ers close to the energy threshold of the telescope, which is ZA dependent. 

425 For some telescope pointing directions the total reconstructed integrated light 

426 of shower images can be reduced by up to ~ 20 %. This is not only the case for 

427 low energies but also for TeV 7-rays. Consequently, if GF effects are not taken 

428 into account the energy of 7-ray candidates from observational data will be 

429 systematically underestimated whereas the 7 efficiency will be overestimated. 

430 Both effects degrade the determination of the flux from a 7-ray source if they 

431 are not properly taken into account in the MC simulation. 



432 7 Conclusions 

433 The results from the MC studies show that the GF can significantly affect 

434 both the shape and the orientation of shower images recorded with an lACT 

435 like MAGIC. Therefore, the orientation discrimination of 7-rays against un- 

436 wanted (hadronic) background can be significantly degraded. It was demon- 

437 strated that the de-rotation of the shower images does not help to recover the 

438 pointing entirely. At most 10 % of the events can be recovered by de-rotation 

439 requiring the knowledge of the impact parameter and energy of the 7-rays. 

440 The influence of the GF also degrades the DISP-estimated arrival direction 

441 of MC-generated 7-rays. Due to the influence of the GF on the development 

442 of EAS the DISP distribution can be significantly elongated perpendicular to 

443 the projection of the GF in the camera. The quality of a sky map is degraded 

444 in a way that a point-like source appears extended unless it is compared to 

445 a proper MC simulation taking into account the trajectory of the source in 

446 the sky. However, the peak of the DISP distribution is always centred at the 

447 source position. 

448 It was also shown that the influence of the GF on EAS can signiflcantly affect 

449 the energy reconstruction and the trigger efficiency for 7-rays. If this effect is 

450 not taken into account, the energy of 7-ray candidates from observational data 

451 will be systematically underestimated (up to ~ 20 % effect). For low energies 

452 close to the analysis threshold (< 100 GeV) the 7 efficiency also depends on 

453 the position of the telescope in the Cherenkov light pool [1]. At higher energies 

454 (~ 300 GeV- 1 TeV), the 7 efficiency is affected only at large ZA, where the 
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455 telescope threshold energy is significantly increased (< 50% effect). 

456 It was demonstrated that the extent of the GF effects not only depends on 

457 the orientation of EAS with respect to the direction of the GF but also on 

458 the position of the telescope with respect to the EAS core location on ground. 

459 Shower images are not only rotated away from the projected direction of the 

460 GF in the telescope camera plane but can also be rotated towards it, contrary 

461 to what was reported in [8j. 

462 Altogether, GF effects on EAS affect the 7-ray sensitivity of an lACT and the 

463 determination of the flux from a VHE 7-ray source. Distinct MC data covering 

464 the same ZA and Az angle range as the observational data being analysed are 

465 required to account for GF effects. 

466 It is remarkable that the GF effects not only occur at very low energies but also 

467 at high energies around 1 TeV. The GF effects are rather pronounced at 7-ray 

468 energies around 450 GeV. The reason for GF effects to occur at high energies 

469 is presumably linked to a characteristic feature in the development of a 7-ray 

470 induced EAS. The process of multiplication in EAS continues until the aver- 

471 age energy of the shower particles is insufficient to further produce secondary 

472 particles in subsequent collisions. At this stage of the shower development, 

473 the shower maximum is reached (largest number of secondary particles) and 

474 the average energy of the secondaries is close to the so-called critical energy 

475 of ~ 100 Me V [26] below which secondary electrons and positrons lose their 

476 energy predominantly through ionisation of air molecules [27]. At the shower 

477 maximum, the average energy of the secondary particles is independent of the 

478 primary 7-ray energy and the GF has on average the same influence on the 

479 secondary particles. 

480 Apart from that, the average atmospheric depth at which the shower maximum 

481 occurs increases logarithmically with increasing energy of the primary 7-ray 

482 [27], and therefore the track along which secondary electrons and positrons 

483 suffer from Lorentz deflection increases, too. 

484 Another point worthy of mentioning is the fact that the threshold energy for 

485 a charged particle to emit Cherenkov light decreases with increasing atmo- 

486 spheric depth. Hence, in high-energy EAS, even charged secondaries of lower 

487 energy suffering strong Lorentz deflection may additionally contribute to the 

488 Cherenkov light pool on ground. 

489 GF effects on the hadron induced background were not studied. It is impos- 

490 sible to show the rotation effect using shower images from hadron candidates 

491 of observational data, because they do not point to any source. Also, possible 

492 GF effects on the hadron induced background presumably do not degrade the 

493 background discrimination. Close to the energy threshold the trigger efficiency 

494 for the hadronic background should be reduced, too. 
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495 8 Outlook 



496 Although this study focuses on GF effects relevant for a single lACT like 

497 MAGIC the influence of the GF on the development of EAS is expected to de- 

498 grade also the performance of stereoscopic lACT arrays. In a stereoscopic tele- 

499 scope system multiple telescopes view the same EAS. The individual images 

500 are then combined to form a common event. Therefore, stereoscopic systems 

501 allow for a three-dimensional reconstruction of the shower axis resulting in an 

502 improved sensitivity. The shower direction on the sky is estimated from the 

503 intersection point of the major image axes in a composite field of view, on an 

504 event by event basis [22] . Due to the rotation of the individual shower images 

505 the intersection point will be modified in a way that the source direction in 

506 the sky is wrongly reconstructed. Stereoscopic systems do not trigger homoge- 

507 neously but preferably on EAS with impact positions between the telescopes. 

508 As the impact distance on ground between the EAS and the individual tele- 

509 scopes of an array is different in the majority of cases, the GF will deteriorate 

510 the orientation of the individual shower images differently. As a result the per- 

511 formance of stereoscopic lACT arrays is degraded by GF effects. A detailed 

512 MC study on the influence of GF effects on the performance of stereoscopic 

513 lACT systems is in preparation. 

514 The intensity of the GF on the Earth' surface ranges from about 20 /iT to 

515 about 70 /xT [T7]. It is therefore important to take into account the GF effects 

516 for the site selection of future projects utilising the imaging air Cherenkov 

517 technique. To minimise the influence of the GF on the detector performance 

518 it is mandatory to select a site with a low absolute value of the GF. Hence, the 

519 best-suited location would be close to the so-called South Atlantic Anomaly, 

520 where the GF strength is minimal, amounting to about one half of the value 

521 for the MAGIC telescope site. 

522 It is difficult to study GF effects in observational data. The elevation effect on 

523 shower images complicates such studies. Directions with strong magnetic field 

524 correspond to large ZA and the sensitivity of an lACT changes as a function of 

525 the ZA as a result of changing shower image characteristics due to increasing 

526 air mass with increasing ZA [22]. There are several requirements a 7-ray source 

527 should fulfil to be an appropriate candidate for GF studies in observational 

528 data: it should be strong, preferably point-like, stable and it should follow a 

529 trajectory corresponding to a large GF component normal to the telescope 

530 pointing direction (figure [1] (a)). 

531 Preliminary results from studies on GF effects in observational data taken with 

532 MAGIC were already shown in [T|28] . It was demonstrated that the pointing 

533 resolution of MAGIC allows to study GF effects in observational data even for 

534 a very low component of the GF normal to the shower direction. However, an 

535 extensive study on GF effects in observational data is in progress. 
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